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Abstract 

A sensitivity study for aircraft wake vortex 
transport has been conducted using a validated large 
eddy simulation (LES) model. The study assumes 
neutrally stratified and nonturbulent environments 
and includes the consequences of the ground. The 
numerical results show that the nondimensional 
lateral transport is primarily influenced by the 
magnitude of the ambient crosswind and is 
insensitive to aircraft type. In most of the 
simulations, the ground effect extends the lateral 
position of the downwind vortex about one initial 
vortex spacing (b„) in the downstream direction. 
Further extension by as much as one b„ occurs 
when the downwind vortex remains “in ground 
effect” (IGE) for relatively long periods of time. 
Results also show that a layer-averaged ambient 
wind velocity can be used to bound the time for 
lateral transport of wake vortices to insure safe 
operations on a parallel runway. 

Nomenclature 

b 0 initial vortex separation 

L lateral distance traveled by the downwind 

vortex 

T nondimensional time; tV o /b 0 

t time 

At time step 

U crosswind component of velocity 

U ma .x ambient crosswind velocity at elevations 

hop and above 

U e layer-averaged ambient crosswind velocity 
U e ambient crosswind velocity 

V 0 initial vortex descent velocity; T 0 / (2 n b n ) 
y fl initial vortex lateral position 

\\irifj drift model lateral position 

y t as s TASS vortex lateral position 

x , y, z longitudinal. lateral, vertical space 
coordinate 

Zi altitude of initial vortex 


Zu>p altitude above which the ambient wind 

velocity becomes constant 
Zwv altitude of vortex core 

a power function in ambient wind profiles 

T 0 initial vortex circulation 

Introduction 

The National Aeronautics and Space 
Administration is making efforts through the 
Terminal Area Productivity Program to safely 
increase terminal capacity and reduce delays to 
meet the increasing air travel demands. 

Particularly, the Aircraft VOrtex Spacing System 
(A VOSS) integrates systems that assess the present 
and short-term weather state and make transport 
and decay predictions for aircraft wake vortices [1, 
2, 3]. These predictions aim to safely decrease 
spacing between arriving and departing aircraft 
under instrument meteorological conditions (IMC) 
and effectively increase airport capacity. The initial 
goal of AVOSS is to successfully demonstrate a 
capability to increase throughput on a single 
runway. Currently, capacity is also reduced during 
IMC when operations on a parallel runway are 
either shutdown or reduced due to the uncertainties 
in the positions of laterally shifting wakes. 

Presently, independent approaches with lateral 
separations of less that 2500 ft (-760 m) are not 
allowed during IMC. This restriction is due to the 
danger of potential encounters with drifting 
vortices. Recent research has been directed at 
improving this limitation without jeopardizing 
safety [4, 5. 6. 7], Schilling [4] concluded that a 
minimum crosswind magnitude of 2.5 m s l was 
necessary to transport vortices to the parallel 
runway (-518 m separation) at the Frankfurt 
International Airport (FIA). This conclusion was 
partly based on numerical simulations and assumed 
a wake vortex lifetime of 180 \. Kopp [5] used a 
continuous-wave lidar at FIA and found that a 
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crosswind of 4 m s J would transport a wake into the 
parallel runway safety area in less than 90 .v. Rudis 
et aL [6] have observed wakes traveling beyond 500 
m. in the presence of strong crosswinds and low 
turbulence. A model that incorporates factors such 
as mean crosswind, runway separation, and wake 
vortex lifetime is needed to insure safe operations 
on parallel runways. 

A wake vortex is considered “in ground 
effect" (IGE) when its altitude is within one initial 
vortex separation ( b 0 ) above the ground. According 
to Critchley and Foot [8], most aircraft wake 
encounters occur between 30-60 m (100-200 ft), 
which for many aircraft is IGE. Wake vortex 
behavior IGE has been investigated both 
observationally [5. 6, 9. 10. 11. 12] and numerically 
[4, 7, 13, 14, 15, 16, 17], but most of the research to 
date has only qualitatively described the 
complexities of this phenomenon. From recent 
large eddy simulation (LES) results, Proctor et ai 
[17] has proposed an empirical model quantifying 
the lateral spreading of a vortex pair once the vortex 
has reached its minimum height in ground effect. 
Their results were limited to cases with little or no 
ambient wind shear. 

Vertical gradients in the ambient crosswind 
have been shown to affect wake vortex descent and 
lateral transport [18, 19, 20]. Specifically, 

nonlinear shear in the crosswind can reduce the 
descent rate, stall, and potentially lead to the rising 
of a wake vortex. The orientation of the wind shear 
gradient (i.e, efU/dz ' positive or negative) will 
determine whether the upwind or the downwind 
vortex is more affected [19]. The combination of 
wind shear and ground effect, which often coincide, 
complicate the problem of predicting wake vortex 
transport. 

The present study will investigate wake vortex 
transport in proximity of the ground in order to 
determine safe bounds for arriving and departing 
aircraft on a downwind parallel runway. The 
simulations conducted in this study assume 
environments with no ambient turbulence. It should 
be noted that turbulence influences vortex decay 
and can have an affect on transport [21]. For 
example, Han et al [22] have examined the 
influence of convective boundary layer turbulence 
eddies on wake vortex transport. They determined 
that substantial lateral shifting of the vortex can 


occur under the influence of turbulence eddies, 
potentially masking the contributions of crosswind 
shear and ground effect. In the present study, 
turbulence effects will be ignored in order to isolate 
the ground effect but should be considered for 
future study. 

In this study we subject a range of aircraft to 
varying ambient crosswind conditions. We will use 
these results to quantify the extent of IGE and 
present a bound on vortex transport to ensure safe 
operations on a downwind parallel runway. 

The Numerical Model and Initial Conditions 

The model employed in this study, the 
Terminal Area Simulation System (TASS) [23], is a 
compressible. LES model, which has successfully 
demonstrated wake vortex behavior in various 
meteorological conditions [16. 17, 19, 22. 24. 25. 

26.27.28.29.30] . 

Grid Configuration 

For the purpose of examining wake transport, 
and to save on computational time, the model is 
integrated in time using a two-dimensional domain. 
The two-dimensional domain inhibits vortex decay 
and prevents three-dimensional instabilities, such as 
Crow linking, from developing [31]. Given these 
limitations, previous studies have shown that 
simulations with the two-dimensional assumption 
have given valid results for vortex transport [16, 19. 

24.30] . 

The coordinate system is chosen such that x is 
along the flight path or longitudinal direction, y is in 
the cross track or lateral direction and z is height 
above ground. The experiments assume a uniform 
one meter grid size. The domain depth is designed 
so that two vortex separations (2 h n ) exist between 
the initial vortex altitude Or,) and the top of the 
domain. The numerical grid translates with the 
downwind wake vortex; therefore, allowing the 
domain width to be independent of crosswind 
magnitude. In all experiments, the width of the 
domain is fixed at 1 1 .5 /?„. 
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Figure 1: Ambient crosswind profiles for TASS 
simulations. 


Atmospheric Conditions 

The simulations are initialized with a neutrally 
stratified and dry atmosphere. The ambient 
crosswind is assumed to be horizontally uniform 
with its vertical distribution prescribed by the 
following power law: 


U,(z)=U~ 



''top 


) 


for k, < Ztop 


for z > Ztop 


( 1 ) 


Three wind profiles with contrasting intensities 
of shear are chosen for the sensitivity study (see 
Fig. 1). Four wind magnitudes are used in each 
profile and are scaled by the initial vortex descent 
rate of the aircraft tested, i.e. U max fV 0 = 1, 2. 4. and 
8. For a typical Boeing 727. these values 
correspond to a dimensional crosswind velocity of 
U,»nx =1-53. 3.06, 6.12, and 12.24 m s' 1 , respectively 
[1 m s 1 = 2 kts]. The parameters for the wind 
profiles are listed in Table 1. 


Table 1. Values for wind profile parameters. 


Profile 

a 

Ztop 

A 

0.4 

10 bo 

B 

0.1 

\0b n 

C 

0.1 

b„ 




Nondimensional Time (t V G / b 0 ) 


Figure 2: Sensitivity to aircraft type for TASS 
simulated trajectories for (a) calm wind and (b) 
crosswind cases. Figure 2b is with wind profile 
“A”, Umax = 2V„ and z, = l.5b 0 - 


Vortex Initialization 

The initial vortex profile is representative of a 
post rollup, wake vortex velocity field. The vortex 
initialization is the same as that described in Proctor 
et ai [17]. The initial vortex parameters for the 
seven aircraft are listed in Table 2 and represent the 
three aircraft weight classes, i.e. small, large, and 
heavy. Six vortex initialization heights are tested: 
Zi/ho = 0.5, 1, 1.5, 2, 4. and 8. The six cases were 
chosen to represent vortices initialized in, near, and 
out of ground effect [32]. 
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Table 2. Initial vortex parameters. 


Aircraft 

b„ ( mi ) 

V 0 (ms 1 ) 

r„ (///V) 

in 

26 

1 .53 

250 

737 

22 

1.49 

205 

747 

50 

1.70 

534 

757 

30 

1.63 

307 

F28 

20 

1.28 

161 

EA 330 

47 

1.28 

378 

DC 10-30 

40 

1.96 

493 


Results 


Sensitivity to Different Aircraft 

To determine if wake vortex transport is 
sensitive to the generating aircraft, two sets of 
experiments are conducted using the seven aircraft 
listed in Table 2. Only the downwind (downshear) 
vortex is analyzed since it is usually the first to 
impact parallel runway operations. 

The first experiment examines the lateral 
transport of the downwind vortex in an atmosphere 
without crosswind (calm wind). A comparison of 
the simulations using conventional nondimensional- 
izations [33] is shown in Fig. 2a. During the initial 
period of the simulation (T < 1) the vortices 
descend due to their mutual induction and exhibit 
no lateral displacement. By T = 2, the vortices lose 
their vertical motion and spread laterally due to 
ground effect [34]. Between T - 3 and 5, viscous 
interaction with the ground causes the development 
of a shear layer with opposite sign vorticity to that 
of the primary vortices. This shear layer separates 
from the surface and is transported around the 
periphery of each of the vortices [9]. Secondary 
vortices are generated within this shear layer and 
orbit the primary vortices. The circulation from the 
secondary vortices causes the primary vortices to 
rise and bounce. Later, the primary vortices begin 
to oscillate under the influence of the orbiting 
secondary vortices. It should be noted that the 2-D 
simulations might exaggerate these oscillations 
[16]. However, comparisons with observed vortex 
behavior have shown good predictions by the 2-D 
simulations [16, 24. 30]. 

The second experiment examines lateral 
transport in an atmosphere with crosswind (Fig. 
2b). Similar to the calm wind cases, the curves for 
lateral transport collapse when nondimensionalized. 


revealing the insensitivity to aircraft type. The 
effects of the ground and secondary vortices on 
vortex motion (Fig. 2a) appear to be less significant 
due to the dominating influence of advection, i.e. 
the process of transport by the crosswind velocity 
field. In fact, by T = 8. the lateral position of the 
vortices is an order of magnitude greater than of 
those in calm wind, thus revealing the dominance of 
advection. 

A “Drift Model” Analysis 

An additional 72 simulations have been 
conducted examining wake vortex transport near 
the ground in a variety of crosswind conditions. 
The simulations involve a mix of the three wind 
profiles (Fig, 1), four wind magnitudes (t/ mAV ), and 
six initial vortex altitudes (c,). which were defined 
earlier. A simple relationship for vortex advection, 
i.e. a drift model, has been calculated for each 
simulation: 

>v(h= y.tnfM - &)+u ,(z m .w (2) 

where U e is from Eq. (1). In this calculation, ~ UT . 
the height of the vortex at time t. is taken from the 
TASS simulation. 

In order to evaluate the effects of ground and 
the nonlinear interaction with a non-uniform 
crosswind profile, the difference is taken between 
Eq. 2 and the lateral position predicted by the TASS 
simulation. The resulting ‘drift differences’ show 
that in nearly all cases. TASS predicts a lateral 
position downwind of the position calculated by the 
drift model (Fig. 3). This downwind extension 
indicates that transport of the downwind vortex IGE 
is greater than that due to advection alone. The 
majority of the ‘drift difference’ curves asymptote 
close to one b () (see Fig. 3). Thus, the ground effect 
extends the lateral position of the downwind vortex 
to about one initial vortex spacing. 

A few cases exhibit a ‘drift difference' 
(extension) greater than one b„. This further 
extension is attributed to either multiple IGE 
encounters or a prolonged ground interaction. 
Figure 4 shows a vortex path revealing a prolonged 
ground interaction. The vortex initially sinks and 
quickly rebounds to a height of 1 .5 b ( >. then it begins 
to oscillate near z = b n . Lidar measurements from 
the Memphis-95 field study provided observational 
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Drift Difference [(y - y )/b ] Drift Difference [(y - y )/b ] Drift Difference [<y - y )/b ] 



N ondim ensional Time (t V G / b„) 

Figure 3: Drift differences between TASS prediction and drift model (Eq. 2) for cases with (a) wind 
profile “A”, (b) wind profile “B”, and (c) wind profile “C”. 








evidence of similar wake paths [35]. Of 80 IGF 
cases identified, 8 (10%) displayed a prolonged 
ground interaction that is similar to that depicted in 
Fig. 4. 

It appears that the prolonged ground 
interaction can be related to the ambient wind shear 
and the initial vortex altitude. For simulations with 
relatively small magnitudes of ambient crosswind 
shear gradient, i.e. \c 2 U/dz 2 \ < 0.001 m l s \ there 
was little affect on the vortex descent rate. In these 
cases, a prolonged ground interaction occurred and 
drift differences were greater than one b n . On the 
other hand, for environments with relatively large 
magnitudes of the crosswind shear gradient, i.e. 
| cU/c^l > 0.001 m l s'\ the descent rate was 
reduced. For these cases the vortices ascended after 
ground interaction, and their drift differences 
remained bounded by one b f> . The vortex 
initialization height can be as important as the 
magnitude of shear. When the vortex was 
initialized IGF, i.e. n = 0.5 h () , a prolonged ground 
interaction occurred regardless of the magnitude of 
the ambient wind shear gradient. In these 
simulations, the vortex eventually ascended above 
z = b a , and the drift difference exceeded one b () . 

The drift model analysis shows that the ground 
effect extends the lateral position of the downwind 
vortex to about one initial vortex spacing (b„), 
relative to transport without ground effect. Further 
extension beyond one b 0 occurs when the 
downwind vortex either remains in ground effect 
for long periods of time or rebounds in ground 
effect more than once. 



Figure 4: Vortex trajectory for the TASS case 
with wind profile “A”, Umax = 4V 0 , and Zi = lb (l . 


A “Time-to-Cross” Analysis 

A drawback from implementing the results of 
the previous section into an operational system is 
knowledge of the vortex altitude. A layer-averaged 
approach is defined and evaluated for a prediction 
of vortex lateral transport, rather than estimating the 
height of the vortex. This approach considers 
neither ground effect nor the vertical trajectory of 
the vortex, and is used to predict the transport of the 
downwind vortex. 


The time interval for the downwind vortex to 
travel across a selected distance is compared 
between the layer-averaged prediction and the 
TASS simulation (Fig. 5). The layer-average 

prediction is determined by t-L /U , where U is 
the layer-averaged ambient crosswind velocity 
(averaged between z - 0.5 h () and z = Zi+ 0.5 b n ). 


U = 


rz t + 0.5 h (i 

J 0,5b o 

[Zj +0.5h t 
J0.5h o 



(3) 


and where L is an arbitrary distance from the initial 
vortex position. The predictions are determined for 
values of L that range from 100 to 800 m by 100 m 
increments. 

An empirical relationship for bounding the 
time interval for lateral transport of the downwind 
vortex IGF is determined from the data in Fig. 5, 
and is: 


/ = 0.688 U\J (4) 

The above relationship (i.e. transport bound. Fig. 5) 
represents the minimum time interval for the vortex 
to transport across distance L (L can represent the 
separation between parallel runways). The thin line 

in Fig. 5 represents t-L /U . which is the transport 
by simple advection. 

This analysis is applied to nine observed wake 
vortices and is shown also in Fig. 5. The observed 
wake positions were measured by lidar during the 
Memphis-95 field experiment [35]. The wind 
profiles used in this analysis were from 
measurements near the lidar site. The predicted 
values for the observed vortices are within or just 
outside of the transport bound. 
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Simulated or observed time (s) 



Figure 5: Time to travel distance L from 
prediction with Eq. 3 versus TASS simulation 
(small symbols). Also included is the prediction 
versus observed time for several Memphis cases 
(large symbols). 

Summary and Conclusions 

Wake vortex transport in the proximity of the 
ground has been examined. The simulations 
conducted assume 1) no ambient turbulence. 2) 
little to no vortex decay, 3) mean crosswind 
profiles, and 4) neutral stratification. Results of the 
simulations are used to evaluate lateral transport 
bounds in ground effect. The conclusions of the 
present study are: 

• Lateral transport is primarily influenced 
by the magnitude of the ambient 
crosswind. 

• Nondimensional lateral transport is 
insensitive to aircraft type. 

• The ground effect extends the lateral 
position of the downwind vortex to about 
one initial vortex spacing (/?„), relative to 
transport without ground effect. 

• Further extension beyond one b <y occurs 
when the downwind vortex experiences a 
prolonged ground interaction. 

• Results of all the simulations show that a 
layer-averaged ambient wind velocity 
can be used to bound the predicted time 
to travel across a specified distance, i.e. t 

= 0.688 L/U . This distance can 
represent the separation between parallel 
runways. 


Future studies should address the sensitivity of 
wake transport to ambient wind shear both in and 
out of ground effect. Also, a select number of cases 
should be extended to three-dimensions to 
determine the effects of decay on wake transport. 
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